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[57] ABSTRACT 

A low-bit-rate, low-cost, all-digital preambleless de- 
modulator for maritime and mobile data communica- 
tions operates under severe high noise conditions, fast 
Doppler frequency shifts, large frequency offsets, and 
multipath fading. Sophisticated algorithms, including 
an FFT-based burst acquisition system, a cycle-slip 
resistant carrier phase tracker, an innovative Doppler 
tracker, and a fast acquisition symbol synchronizer, 
provide reliable burst reception. The compact DSP- 
based demodulator includes an input buffer receiving a 
complex sampled baseband input signal and providing a 
baseband output to a coarse frequency estimator fast 
Fourier transform (FFT) or discrete Fourier transform 
(DFT) module which produces a first estimation of the 
carrier frequency. A fine frequency estimator FFT or 
DFT module receives the first estimation and provides 
a second estimation of the carrier frequency. An extra 
coarse frequency estimator FFT or DFT module may 
be provided between the buffer and the coarse fre- 
quency estimator. 

18 Claims, 11 Drawing Sheets 
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This is a continuation-in-part of application Ser. No. 
07/800,020, filed Nov. 29, 1991 and now abandoned. 

BACKGROUND OF THE INVENTION 

The present invention relates to a digitally imple- 
mented demodulator for maritime and mobile data com- 
munications, and in particular to a fast frequency esti- 
mator which operates without using an acquisition pre- 15 
amble. 

A typical packet format used for burst communica- 
tions and signaling in a time division multiple access 
(TDM A) network employs a structure that includes an 
acquisition preamble at the start of the packet, followed 20 
by a known unique word (UW) pattern. The data por- 
tion of the packet then follows, with additional framing 
bits inserted periodically for long packets, with the 
packet ending in another known "end-of-packet" se 
quence. 

The acquisition preamble that precedes the data por- 
tion of the packet typically incorporates an unmodu- 
lated carrier sequence for the carrier frequency and 
phase estimate at the receiver, followed by a clock 3Q 
recovery sequence for proper receiver clock phase 
alignment. The UW pattern is used for phase ambiguity 
resolution and burst time synchronization. After the 
preamble and UW segments have been received and 
acted upon, the receiver is ready to demodulate the 35 
ensuing data segment with the correct frequency, 
phase, and clock adjustments. 

The unmodulated carrier segment appears as a single 
tone in the frequency spectrum for a short duration. 
The detection of the carrier can therefore be accom- 40 
plished by frequency domain acquisition algorithms. 

In an analog implementation, a frequency lock loop 
(FLL) may be used. A bank of analog bandpass filters 
followed by energy detectors may be used to give a 
coarse estimate of carrier frequency. A similar tech- 45 
nique utilizes a bank of correlators tuned to several 
discrete carrier frequencies. In a digitally implemented 
receiver, a DFT/FFT-based algorithm may be used. 
Time domain techniques may also be used for the car- 
rier frequency estimate, whereby the unmodulated car- 50 
rier phase differences are computed periodically and an 
average estimate is obtained that is indicative of the rate 
of the phase change or the frequency offset. 

The methods used for carrier frequency acquisition 
described above have limitations that preclude their use 53 
for communication networks that operate under hostile 
channel environments that include low signal-to-noise 
levels, large frequency offsets, Doppler frequency 
shifts, and multipath fading. The FLL technique, in 
particular, requires a finite "lock-in" time for the loop to 60 
acquire. The tracking range is also very limited if a 
high-resolution measurement is required. For data pack- 
ets with very short acquisition preambles, or for pream- 
bleless packets, the FLL technique is disadvantageous. 
The FLL method may still be used during the data 65 
portion of the packet if the signal is squared to remove 
the modulation (for binary phase shift keying, or BPSK) 
but with a corresponding loss of operating signal-to- 



noise ratio that will compromise its operation for low 
channel signal-to-noise ratios. 

The bandpass filter and correlator approaches be- 
come excessively complex if a very high resolution 
measurement is required. If the filters or correlators are 
of the analog type, they are also subject to drift and 
require precise components and frequent adjustments. 
In addition, signal detection is sensitive to the input 
signal level unless an automatic gain control (AGC) 
circuit is employed. 

The straightforward FFT technique suffers the disad- 
vantage of susceptibility to multipath that will affect the 
signal strengths in the frequency bins. A single DFT or 
FFT must also be very large (1024 points or greater) to 
operate successfully at very low values of Es/ No. A 
threshold to detect the signal presence will be unreliable 
unless external AGC is employed that will track the 
multipath fading. A tradeoff DFT/FFT computation 
time versus frequency resolution is also necessary. Fine 
frequency resolution translates to longer FFT computa- 
tion time, and vice" versa. In addition, this technique is 
suited for the case of a packet with an acquisition pre- 
amble. For preambleless packets, the straightforward 
FFT technique requires modifications. 

The time domain frequency estimation technique is 
dependent on the unmodulated portion of the acquisi- 
tion preamble, and is totally unsuited for preambleless 
packets. This method is also particularly sensitive to the 
magnitude of the frequency offset. If the frequency 
offset is large enough to cause a rotation of the signal 
vector through 2tt radians between the periodic esti- 
mates, then the measurement will be in error. In addi- 
tion, due to the low operating signal-to-noise ratios, a 
long time average of the phase difference measurements 
are required, thus necessitating a long acquisition pre- 
amble. 

The required preamble that precedes the data portion 
of the packet often constitutes an excessive overhead 
that reduces the channel transmission efficiencies for 
short data packets. Particularly in the case of mobile 
and maritime communication networks that experience 
signal fades of varying degrees, the system design man- 
dates the use of very short data bursts (approximately a 
few hundred bits) that are relatively immune to fades. 
These bursts generally are used for the channel request 
and assignment functions from the remote terminal to a 
central location, typically a network coordination cen- 
ter. With additional channel impairments such as Dop- 
pler shifts, low carrier-to-noise ratios, and frequency 
offsets, long acquisition preambles usually are required 
to locate the carrier correctly, and to acquire its fre- 
quency and phase. The use of a long preamble in this 
case is an extremely undesirable overhead that seriously 
undermines the channel transmission efficiencies for a 
large communication network, with potentially thou- 
sands of such remote terminals. 

SUMMARY OF THE INVENTION 

In view of the foregoing, it is an object of the present 
invention to provide a compact, low bit rate, low cost, 
all-digital burst demodulator that requires absolutely no 
acquisition preamble. Both the acquisition algorithms 
and the hardware architecture are applicable under a 
wide variety of system constraints, and are particularly 
well suited for reception of short data packets under 
very adverse channel conditions. 

In accordance with the present invention, a digitally 
implemented fast frequency estimator has been devised 
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for application in low bit rate maritime and mobile 
MPSK-modulated (M'ary phase shift keying) data com- 
munication that performs reliably under severe channel 
conditions, such as low signal-to-noise ratios, large fre- 
quency offsets and multipath fading. The estimator is 5 
well suited for carrier frequency estimation for short 
burst packets, particularly for those packets that do not 
have an acquisition preamble preceding the data por- 
tion, such as the preambleless packet used for the IN- 
MARSAT Standard-C signalling channel. The estima- 10 
tor detects signal presence independent of the input 
signal level, eliminating the need for automatic gain 
control. The estimation technique is particularly suited 
for implementation in a digital signal processor. The 
limitation on the operating bit rate is therefore solely a 15 
function of the processor speed. 

In implementation, a two step process, using Discrete 
Fourier Transforms (DFT) or Fast Fourier Transforms 
(FFT), is used. The first step involves taking a time 
average of a sequence of hopping FFTs, with the sam- 20 
pie window slid sequentially by half the window length 
for each such DFT/FFT computation. A time average 
of the signal energy in each of the DFT/FFT frequency 
"bins", followed by a relative energy level comparison 
between the bins, yields (a) a detection of the burst 25 
presence, and (b) the location of the carrier and its offset 
frequency from the reference frequency. The coarse 
carrier frequency is resolved to ±75 Hz in the first step. 
The second step involves performing another 
DFT/FFT around the frequency range identified by 30 
the first step to yield a fine frequency estimate resolved 
to less than 10 Hz. 

The present invention embodies several improve- 
ments and enhancements over the old methods de- 
scribed above. A summary of the advantages is as fol- 35 
lows: 

1. The estimator operates on MPSK modulated data 
directly, eliminating the requirement of a carrier 
acquisition preamble sequence; 

2. The use of separate coarse and fine estimators al- 40 
lows carrier acquisition over a very wide fre- 
quency range (greater than the MPSK symbol rate) 

at very high resolution (< 1% of the symbol rate). 
This resolution accuracy, under conditions of large 
frequency offsets, low signal-to-noise ratios and 45 
multi path f ading, is a major improvement in the 
DFT/FFT-based method. In one embodiment of 
the invention, an additional estimator stage is 
added, yielding a multi-stage approach which ena- 
bles yet further improvement; 50 

3. Time averaging of the measured signal strengths in 
each of the DFT/FFT frequency bins increases the 
signal-to-noise ratio of the measurement, making its 
operation possible even at very low channel signal- 
to-noisc values; 55 

4. The use of "time-hopped" DFT/FFTs with the 
sample window moved by half the window width 
for each DFT/FFT computation provides the ad- 
vantage of minimizing the variance of the estima- 
tion while drastically reducing the number of com- 60 
putations as compared to DFT/FFTs slid in time 
by 1 sample period; 

5. Relative energy level comparison between 
DFT/FFT bins removes the effect of signal 
strength variations due to multipath fading and 65 
input gain variations. Signal activity is detected 
independent of the input signal level. As a result, 
there is no need for an AGC as otherwise typically 



would be employed. The result of performing op- 
erations (3), (4), and (5) yields not only the detec- 
tion of the burst presence, but also locates the car- 
rier and determines its offset frequency from the 
reference frequency; 
6. The use of sample rate decimation between stages 
greatly reduces the computational complexity 
compared to other FFT based approaches. Actual 
implementation of this algorithm has been opti- 
mized for digital signal processor-based (DSP) 
demodulators. Thus, with the DSP assembly lan- 
guage code already developed, an increase in the 
operating bit rate is readily achieved by migrating 
the code to faster processor versions in the same 
DSP family or to another faster processor of a 
different family. 
The use of the second DFT/FFT for the fine fre- 
quency estimation and the application of this technique 
for the coarse and fine frequency estimations for pream- 
bleless data packets are particularly important features 
of the inventive method. The use of multiple DFT/FFT 
stages for carrier frequency estimation, in which the 
output of each preceding stage is fed to the next stage, 
greatly increases the frequency resolution at each stage. 

U.S. Pat. No. 4,466,108, commonly assigned with the 
present application, discloses a technique for perform- 
ing phase shift keying (PSK) synchronization without 
requiring a preamble. However, the technique disclosed 
therein does not use the multiple-stage DFT/FFT ap- 
proach of the present invention, inter alia, as will be 
discussed herein. 

Other known techniques, such as that disclosed in 
U.S. Pat. No. 4,245,325, use an FFT technique, but for 
a very different purpose from that of the invention. U.S. 
Pat. No. 4,912,422 uses multiple FFT stages, but uses 
the resulting frequency estimate to steer a separate set of 
bandpass filters, whereas in the present invention FFT 
and filtering processes are combined in a single opera- 
tion. 

Other patents considered to be possibly relevant 
background to the present invention include U.S. Pat. 
Nos. 4,528,567; 4,618,830; 4,654,667; 4,689,806; 
4,827,488; 4,870,420; and 4,885,756. 

Two embodiments of the software algorithms to be 
implemented in the hardware embodiments of the in- 
vention are described herein. The second embodiment 
has greater flexibility than the first in terms of operabil- 
ity at different symbol rates. In particular, the second 
embodiment provides yet another, extra coarse resolu- 
tion DFT/FFT stage, yielding a total of three 
DFT/FFT stages. The addition of the "extra coarse" 
stage estimates the frequency of the modulated spec- 
trum. The input signal is not preprocessed by a squaring 
operation, as has been done traditionally to remove the 
modulation first. As a result, performance at low S/N 
ratios is improved greatly, since squaring can degrade 
the S/N ratio by more than 6 dB. The first, "extra 
coarse" stage acts to filter the input signal to remove a 
significant amount of noise. Squaring then is employed 
in the second stage, where the impact of the degrada- 
tion is not as great. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other features and advantages of 
the invention will be understood from the following 
detailed description taken in conjunction with the ac- 
companying drawings, in which: 
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FIG. 1 shows a preambleless packet format for the testing of the invention was carried out in accordance 

Inmarsat -C signaling channel, at two different symbol with these parameters. 

rales ; These constraints impose very challenging process- 

FIG. 2 is a functional block diagram in accordance ing requirements for the preambleless demodulator in 

with a first embodiment of the present invention; 5 order to successfully receive and detect each preamble- 

FIG. 3 is a block diagram of the Doppler estimator less packet under the worst possible channel environ- 

and frequency tracker of FIG. 2; ment and slot timing uncertainties. As now will be de- 

FIG. 4 is a block diagram of a phase tracker and bit scribed, sophisticated acquisition algorithms that func- 

detector used in connection with the Doppler estimator tion in this hostile channel environment have been de- 

and frequency tracker of FIG. 3; 10 velopcd and implemented successfully in DSP-based 

FIG. 5 is a structural hardware block diagram of a hardware. Test data obtained with the hardware imple- 

first embodiment of the present invention; mentation and a specialized test fixture indicate excel- 

FIG. 6 is a block diagram of a test setup used to test lent correspondence with the expected acquisition per- 

the algorithms used in the preambleless demodulator formance, as predicted by high-level computer simula- 

according to the present invention; 15 tion as well as theory. These results will be summarized 

FIG. 7 is a functional block diagram of the inventive below, 

frequency acquisition algorithm implemented in Inmar- In the following, it is assumed that an MPSK modu- 

sat-C Coast Earth Station (CES) signalling/message lated signal has undergone IF conversion to baseband, 

demodulators in accordance with the first embodiment; and that both the in-phasc and quadrature (I and Q) 

FIG. 8 is a functional block diagram of a frequency 20 components are available to the digitally-implemented 

acquisition algorithm in accordance with a second em- demodulator. The baseband received signal will contain 

bodiment of the present invention; a main spectral lobe within the frequency range ±R 5 , 

FIG. 9 is a functional block diagram of a frequency where K s is the symbol rate. Given a large frequency 

acquisition algorithm in accordance with a modification uncertainty, ±f f , in the received carrier due to Doppler 

of the second embodiment of the present invention; shifts and downconverter/upconverter translation er- 

FIG. 10 is an overall block diagram of a three-stage rors, the total receiver one-sided bandwidth must be 

frequency estimator in accordance with the second greater than (R,*ff c ). This represents a lower bound on 

embodiment of the invention; and the receiver input filter bandwidth. In general, how- 

FIG. 11 is a functional block diagram of a frequency ^ ever, the channel spacing determines the filter band- 
acquisition algorithm in accordance with the embodi- width. The receiver I and Q channel filters must exhibit 
ment of FIG. 10. a very sharp rolloflf in order to minimize adjacent chan- 
riPTAn P n nccroTDTinM ncruc nel interfcrence - The frequency acquisition technique 
J^J^cd ^ S^III?^ OF THE described herein assumes the presence of such channel 
PREFERRED EMBODIMENTS fii ters l0 i imit ^ e O ut-of-oand noise and adjacent chim- 
in order to illustrate the inventive frequency acquisi- ne ^ interference. Since the 1 and Q analog-to-digital 
tion algorithms and their hardware implementation, the (A/D) converter modules preceding the digital demod- 
preambleless short data packet format from a new mari- ulator require anti-aliasing filters, this is not a problem, 
time communication network, designated as "Standard- The channel filters also serve as anti-aliasing filters. 
C" or "Inmarsat-C", which has already been specified Tne fi fSt ste P in tn e frequency acquisition method 
by the International Maritime Satellite Organization, consists of computing a series of time-hopped win- 
(INMARSAT), will be used. The preambleless short dowed DFT/FFTs (called "periodograms"), and aver- 
data packet format is used for the slotted ALOHA aging the energy magnitudes over the duration of a 
random access signalling requests from the remote ship burst to reduce noise and fading effects. Given a se- 
terminals. The proposed frame format for this signalling 45 quence of samples, x(n), it is possible to form K sub- 
packet is shown in FIG. 1. This particular frame format sequences x/n), r= 1,2, . . . , K, of length L spaced L/2 
and the Inmarsat-C channel parameters have been used samples apart. For each sub-sequence, the windowed 
as a representative mobile environment to design and DFT/FFT XXm) is given by the following equation: 
demonstrate the efficacy of the preambleless frequency 
acquisition algorithms that were developed. 50 jL-1 

The selected representative preambleless data packet ** m) = 0 x ^ lt ^~ A2n/L) ' tm 
has the following characteristics, as obtained from the 

Inmarsat-C System Definition Manual, for the second where w(n) is the window function. 

generation format shown in FIG. 1 (having a rate of The periodogram I^f m ) is equal to the squared mag ni- 

1200 symbols/sec): No acquisition preamble of any sort; 55 tude of an DFT/FFT bin: 

the TDM A data packet starts with a 64 bit unique word; 

and the data segment is 252 rate-J convolutionally en- WmblWI 2 

coded bits, to yield a packet size of 3 16 bits (632 for the 

600 symbol/sec case). The timeslot width is 370 symbol where 

periods (740 for the 600 symbol/sec case). The severe 60 

channel operating constraints include BPSEC modula- m 

tion at the low rate of 600 and 1200 symbols/sec; a fm " T ~ «™«iu«cy 

carrier frequency uncertainty of ±1450 Hz; a short 

term Doppler rate of 65 Hz/sec with a maximum range The unnormalized power spectral density is then: 

of ±50 Hz; an operating Ej/N 0 of4.7dB at the demodu- 65 

lator input; slot timing uncertainty of ±27 symbols; and # 

finally, a carricr-to-mullipath ratio of 7 dB with a Rician SAfm) - ^1 ijj m ) 

fading bandwidth of 0.7 Hz. As will be seen below, 



03/01/2004, EAST Version: 1.4.1 



5,272,446 

1 8 

The effect of the modulation is first removed by tions need to be taken in minimizing the truncation error 

squaring the complex inputs (for BPSK modulation). while calculating the argument 2fffnT. 

This concentrates most of the signal energy at a single Looking now at a hardware implementation of the 

frequency equal to twice the modulated carrier fre- invention, as seen in FIG. 2, the DSP-based demodula- 

quency. Although the squaring of the signal results in a 5 tor accepts a complex sampled baseband input either 

loss in the operating signal-to-noise ratio, this loss is from a pair of A/D converters connected to the IF 

regained due to the averaging process performed on the downconverter unit, or from a personal computer capa- 

signal strength estimates. Only the 256 symbols (2048 ble of generating simulated real-time channel packet 

samples) in the middle of a timeslot assumed to contain data for test purposes. 

the carrier are processed. The remaining symbols at the 10 An input buffer 205 stores complex input samples for 

beginning and end of the timeslot are discarded due to an entire burst duration before processing begins, thus 

burst location uncertainty. At a sampling frequency of allowing the burst to be reprocessed during several 

fi-9600 samples/sec, a 64 point DFT/FFT provides a passes to recover the final data. The buffer is controlled 

resolution of 150 Hz. The resolution is increased to 75 by an external slot timing pulse which marks the begin- 

Hz since the estimate of the squared data must be di- 15 ning of a burst timeslot. The actual burst location can be 

b y twa anywhere within the first 54 symbols allowed for guard 

The 2048 data samples are divided into 64 blocks of time. The Inmarsat-C slotted ALOHA protocol pres- 

64 samples each (K = L=64), with each block overlap- ents several difficult problems for burst acquisition, as 

ping adjacent blocks by 32 samples. Each block is multi- was outlined above. The usage of each burst timeslot in 

phed by a 64 point Blackman window. A Blackman 20 the signaling channel is random, and a given timeslot 

window has a very broad transition width so th at a may be empty or may contain a collision from two or 

earner located at the very edge of an DFT/FFT fre- more remote transmitters. 

quency bin will be attenuated only 1.1 dB, as compared The first step in burst acquisition is estimation of the 

to 3.9 dB for a rectangular window. The carrier fre- carrier frequency offset. The frequency acquisition al- 

quency is determined by choosing the bin containing 25 gorithm is the most critical to successful demodulator 

the maximum power spectrum. Because the range of the performance since it must also detect signal presence at 

frequency search is limited to ±1450 Hz, the a very low signal to noise (S/N) ratio (1.5 dB for an 

^rT^TfT al 8° ntnm be educed to only 40 E,/N 0 of 4.7 dB) because of the large extra input band- 

DFT/FFT outputs. The entire calculation is very effi- width required to allow for frequency uncertainty. A 

cient when written with in-line code in a DSP chip 30 simple energy detector is ineffective for burst presence 

because the small DFT/FFT size allows exclusive use detection since the signal and noise energies are nearly 

of the chip's internal data RAM. The total number of equal. Frequency acquisition is performed in two steps 

processor instruction cycles is approximately 175,000. Coarse acquisition in coarse frequency estimator 220 

The frequency estimator will locate a peak regardless estimates the carrier frequency to within 75 Hz, and fine 
of the presence of a carrier. Burst presence can be de- 35 acquisition in fine frequency estimator 230 narrows the 
tccted by computing the peak/average energy ratio, uncertainty down to 10 Hz. The data in the input buffer 
The peak energy is simply the energy of the periodo- 205 is corrected for frequency offset using a complex 
gram containing the maximum energy. The average multiplication by the sine and cosine of the carrier fre- 
energy is the average of the other 39 periodograms. If quency generated from look-up tables, 
the ratio is greater than a predetermined threshold, then 40 Symbol timing adjustments are performed digitally in 
a burst has been detected. This method is very reliable software. A conventional A/D sample clock adjust- 
over a wide dynamic range, eliminating the need for ment circuit cannot be used, as an entire burst is sam- 
gain control. Simulations indicate that a threshold near pled before symbol timing is acquired. The frequency- 
2.0 produces good results down to EyN 0 =2 dB. corrected samples are first match filtered. Symbol tim- 

The accuracy of the estimate is increased to within 45 ing then is acquired, and the match filter output delay 

several Hz (depending on the amount of Doppler shift line is shitted to the proper timing phase and decimated 

across a burst) using the fine frequency estimator. The to the symbol rate of 1200 symbols/sec 

40 complex outputs of all 64 DFT/FFTs are stored in After symbol timing comes the removal of the burst 

^CT^rrT- 64 Va,UCS corres P° ndin g to tne time uncertainty, correction for the Doppler shift, and 

DFT/FFT bin containing the carrier can be viewed as 50 detection of the phase-coherent data. The burst location 

a decimated time series of the original input, with a is acquired with a non-coherent unique word (UW) 

decimation factor of 64. If the carrier is located in bin correlator 250. Coherent correlation is not possible 

M, the time series X,(M) is a bandpass filtered version because at the maximum Doppler rate of 65 Hz/sec the 

of Xr(n) filtered around the carrier frequency by the carrier phase will drift significantly over the duration of 

windowed Founer coefficients. Thus, if the series is 55 the UW. The UW correlator 250 serves a second func- 

processed by another L-point DFT/FFT and the bin tion in addition to slot timing acquisition. If the magni- 

with the maximum squared magnitude is selected, the tude of the correlation exceeds a predetermined thresh- 
resolution of the estimate can be increased by a factor of old, then a valid packet has been received Otherwise a 

L. For the representative case of L=64, the fine fre- packet collision is declared and the data is ignored 

quency estimator resolution is 1.2 Hz, assuming no car- 60 Coherent data detection is not possible unless the 

ner frequency drift over the measurement interval. For Doppler shift is removed. The Doppler estimator/- 

mobile satellite applications, the actual resolution is smoother 260 calculates a frequency offset for each 

limited by short-term Doppler effects and is on the symbol and correction is accomplished with a complex 

order of 5-10 Hz. multiplication stage. The Doppler-corrected symbols 

Demodulation of the stored complex samples by the 65 will have arbitrary carrier phase, and the phase estima- 

estimated frequency can be performed by multiplying tor 270 measures the phase offset and adjusts each sym- 

the samples by [cos(27rfnT)+ jsin(27rfnT)]. The sine and bol to the zero degree phase reference. The estimator 

cosine values arc obiained from a look-up table. Precau- 270 is able to detect approaching cycle slips and com- 
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pensate for them before they occur. A coherent UW 
correlation is performed on the output of the phase 
detector to correct for the inherent 180 degree phase 
ambiguity of BPSK modulation. The final data, along 
with status information, is sent to the host processor, 
where Viterbi decoding, data formatting, and cyclic 
redundancy check (CRC) are performed. Further de- 
tails of the major demodulator algorithms are provided 
below. 

With respect to frequency acquisition, the require- 
ment of detecting signal presence/packet collisions as 
well as accurately measuring the carrier frequency of a 
BPSK-moduIated signal at an input S/N near 0 dB is 
especially severe. The inventive DFT/FFT-based spec- 
tral estimation technique was found to be very robust 
and readily suited for DSP implementation. As men- 
tioned above, the method consists of computing a series 
of time-hopped windowed DFT/FFTs and averaging 
the energy magnitudes over the duration of a burst to 
reduce noise and fading effects. The DFT/FFT bin 
with the maximum average energy is presumed to con- 
tain the carrier. A burst detection algorithm, to be de- 
scribed below, detects the presence of the burst from 
the energies in the DFT/FFT bins. Further processing 
using a similar DFT/FFT technique increases the fre- 
quency resolution to well within the Doppler tracker 
range. The method is very reliable over a wide dynamic 
range. An optimal threshold has been obtained through 
simulation, as will be described, yielding excellent re- 
sults. 

As for matched filtering, the lack of any kind of pulse 
shaping at the transmit side results in a simple integrate 
and dump matched filter structure at the receive side. 
Since the symbol timing has not yet been determined, 
demodulated samples are processed through this filter 
with the outputs computed at the sample rate. Once the 
symbol timing has been established, the filter output is 
decimated to the symbol rate. 

Symbol synchronization is achieved by first calculat- 
ing a timing error vector W between the signal enve- 40 
lope and two suitable selected quadrature symbol tim- 
ing references, P(n) and Q(n). The error vector, 
W = U+jV, is then used to estimate the timing phase 
difference by calculating arctan (V/U). The real and 
imaginary components of the timing error vector, U 45 
and V, are obtained from correlating the envelope of 
the matched filter output with the two timing refer- 
ences. After the symbol timing has been calculated, it is 
rounded to the nearest timing phase step, and the proper 
output of the matched filter is selected. 

The symbol rate stability of the transmitter is given as 
±1 part per 10 5 . Therefore, under worst case condi- 
tions, the symbol rate could change by 10~ 3 degrees 
during one signaling packet interval. This is negligible; 
hence, a single symbol timing estimate is sufficient for 
each packet. 

UW correlation is performed after symbol timing has 
been established, but before removal of Doppler fre- 
quency offset and carrier phase recovery. The 64 data 
symbols to be correlated with the UW pattern are di- 
vided down into several subgroups. The correlation 
results of these subgroups with the proper part of the 
UW sequence arc then added and used to determine the 
location of the UW sequence in the packet. 

The Doppler tracker (FIG. 3) must reduce the carrier 
frequency error to less than 2 Hz to prevent a complete 
phase reversal in 3 16 symbols. The tracker removes any 
residual frequency error from the frequency acquisition 
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stage as well as removing the Doppler. The tracker is 
essentially open-loop, to prevent error propagation 
caused by impulse noise. Because the Doppler rate is 
well behaved, only a small number of estimates are 
computed using contiguous N-symbol blocks. N is suffi- 
ciently large to reduce the variance of the estimate to 
much less than 1 Hz. Although the Doppler may drift 
significantly between two estimates, it is assumed that 
the Doppler shifts undergo sinusoidal changes. There- 
fore, the Doppler frequency associated with each sym- 
bol can be computed by an appropriate smoothing func- 
tion between two estimates. 

Two separate estimators are employed for increased 
accuracy. The complex symbols output from the input 
processing block are applied to both estimators. The 
first, coarse estimator 220 employs a discrete Fourier 
transform (DFT)-based technique to estimate the Dop- 
pler frequency. A second, fine estimator 230 is em- 
ployed to substantially increase resolution. The esti- 
mates are sent to the post -processor block, and a new 
Doppler phase angle is calculated for each symbol. The 
phase updates drive a numerically controlled oscillator 
which is used for the Doppler correction of the 
matched filtered symbols. 

The Doppler post-processor ensures that the Doppler 
corrections do not suffer from appreciable phase jitter. 
The phase tracker (FIG. C) adjusts the Doppler-cor- 
rected symbols to the zero-degree phase reference in 
order to perform bit detection. A very simple algorithm 
is employed that is essentially open-loop, again provid- 
ing very fast recovery from noise bursts. The complex 
symbols are squared to remove modulation effects, and 
then the real and imaginary components are averaged 
over a block of ±M symbols to give an average two-di- 
mensional carrier phase coordinate. The carrier phase 
can be derived by making an arctangent calculation 
using a power series approximation. The input symbols 
are then rotated by an NCO so that one of the two 
BPSK modulated symbols falls at zero degrees. A 0/1 
bit decision is made based on the sign of the real compo- 
nent of the rotated symbols. 

The carrier phase normally rotates quite slowly, but 
occasional corrections must be fed back to the Doppler 
tracker to maintain a phase less than 1 80 degrees, or else 
cycle slips would occur. These corrections ensure that 
the first symbol of the next block will have a phase 
offset near zero degrees. The final bit decision is passed 
to the host processor for subsequent Viterbi decoding. 
— FIG. 5 is a high-level block diagram of the DSP- 
based demodulator hardware. The demodulator oper- 
ates on baseband received digitized inputs obtained 
from the output of the IF down-converter or a personal 
computer test interface. High speed erasable program- 
mable read-only memory (EPROM) program memory 
310 and random access memory (RAM) data memory 
520 are used so that the external memory is addressed 
without wait states. Extensive use of programmable 
logic devices yields a compact design for the DSP pe- 
ripheral logic, such as the address decoders, memory 
control signals, port decoders, and I/O interfacing. 
Such implementation is well known to those of working 
skill in this technological field, so further detailed de- 
scription of this aspect of the invention is not believed 
necessary. 

Communication between the demodulator 500 and 
the host processor 950 takes place through a high-speed 
global RAM 525 connected to each processor bus. Out- 
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put ports for data and other status monitoring and con- 
trol signals are provided. 

To test the inventive demodulator, a special personal 
computer-based test setup has been designed using the 
simulated test data with channel impairments. As shown 5 
in FIG. 6, the test files, generated on a larger, minicom- 
puter or mainframe computer (a VAX-1 1/750 is identi- 
fied in FIG. 6), were downloaded to a personal com- 
puter. User selectable frequency offsets, Doppier shifts, 
noise levels, delays, and symbol timing offsets within 10 
the slot could be entered on the larger computer during 
file creation. Test program modules were written on the 
personal computer to format the data file and download 
the sampled values to the demodulator board via a par- 
allel port interface. Demodulated data then was 15 
uploaded from the DSP emulator to the personal com- 
puter through a serial port, where another program 
compared it with the transmitted packet and reported 
the error statistics. 

The preambleless signaling demodulator was tested at 20 
a symbol rate of 1200 bit/sec for frequency offsets rang- 
ing between ±1450 Hz, sinusoidal Doppier shifts of 65 
Hz/sec with a range of ±50 Hz, E s fN 0 values from 2 to 
6 dB and the reference noiseless case, symbol timing 
offsets between ±180 degrees in 45 degree steps, and 25 
slot timing uncertainty in the range of ±27 symbols. 
With no noise, the inventive demodulator recovered the 
data successfully, with no errors, under all of the above 
conditions. At the specified operating E*/N 0 of 4.7 dB, 
an average of 2.7 raw errors per packet (316 symbols) 30 
were detected. This result compares well with the limit- 
ing theoretical value of 2.4 errors per packet, represent- 
ing performance within 0.2 dB of the theoretical bit 
error rate (BER) curve. The error statistics indicate that 
no packet errors will be experienced after Viterbi de- 35 
coding (to be performed by the host processor). Repre- 
sentative test results are provided in Table 1. 

TABLE 1 
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be seen that the invention also has wider ranging appli- 
cations. 

With future enhancements to DSP technology and 
repackaging for medium and high-volume production, 
the entire preambleless demodulator can be produced 
easily in a very compact package at even higher symbol 
rates. Software changes also can be introduced to ex- 
tend the operating bit rate, change the modulation for- 
mat, or customize the operating parameters for the 
Doppier roultipath fading according to the desired 
channel conditions. This technology also lends itself to 
other low bit rate systems, where the existing preambles 
may now be shortened or discarded altogether, if the 
digital implementation approach described herein is 
adopted. 

The frequency acquisition algorithm just described 
for the Inmarsat-C CES signalling and message channel 
demodulators was optimized for operation at 1200 bit/- 
sec, representing the second generation of service. With 
the first generation service introduction at 600 bit/sec, 
the inventors investigated whether the inventive tech- 
nique would still function reliably at the lower rate. It 
was found that in the 600 bit/sec mode the current 
algorithm would still work reliably at an Es/N 0 of 4.7 
dB, but with a margin of less than 1 dB. Therefore, in 
order to operate reliably at 600 bit/sec in presence of 
multipath fading, an enhancement to the current fre- 
quency estimation algorithm was considered to be de- 
sirable. 

In view of the foregoing, in accordance with another 
embodiment of the present invention, another enhanced 
frequency acquisition algorithm that addresses this issue 
is discussed herein. The reliable operation of the new 
algorithm for the first generation, at an Ef/N 0 of 0 dB 
has been validated. This technique can also be applied 
to the second generation service at 1200 bit/sec for 
improved acquisition performance. 



Es/No 
dB 



Representative Measured Test Results for the Preambleless Demodulator 
Input Parameters Measured Values 



fs(Hi) 



Doppier 
on/off 



AT (deg.) 



Slot Delay 
sym 



fs(Hz) AT (deg.) 



Slot Delay 
sym 



No Noise 0 

No Noise + 1450 

No Noise - 1450 

No Noise +1450 

No Noise -1450 
4.7 0 

4.7 +1450 



Errors per 
packet (in 
symbols) 



off 
off 
off 



4.7 
4.7 
4.7 
3.0 
2.0 



-1450 
+ 953 
+ 1450 
-1450 
+ 725 



0 

+ 180 
+ 180 
+ 180 
+ 180 
+ 160 
+45 
-45 
-180 
0 

-180 
0 



4 
3 
3 
3 
3 
3 
3 
3 
4 
53 
3 

53 



0 

+ 1450.75 
-1450.81 
+ 1460.13 
- 1439.06 
+ 9.38 
+ 1460.13 
+ 1441.44 
+963.25 
+ 1460.13 
-1439.06 
+ 733.56 



0 

+ 180 
+ 180 
+ 180 
+ 180 
+ 180 
+ 45 
-45 
+ 180 
0 

-180 
0 



4 

3 
3 
3 
3 
3 
3 
3 
4 

53 
3 
53 



None 
None 
None 
None 
None 

4 

2 

1 

3 

4 
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As has been described, in accordance with the pres- 
ent invention, a compact, low cost, low complexity 55 
DSP-based preambleless demodulator employs power- 
ful digital signal processing techniques to detect and 
acquire the carrier frequency offset and symbol timing 
under high noise conditions, as well as to track Doppier 
and short term variations in signal amplitude because of 60 
multipath fading. As a result, it is possible to acquire and 
demodulate the short burst packets without any acquisi- 
tion preamble. It is noted that, while the invention has 
been described with reference to mobile and maritime 
satellite and terrestrial communication networks that 65 
use very short burst packets for access request and 
response channels (an application for the maritime In- 
marsat-C system having been described in detail), it can 



The Inmarsat-C CES signalling and message demod- 
ulators operate on baseband complex samples. With 
BPSK modulation and no pulse shaping present on the 
link, the one-sided bandwidth of the input signal spec- 
trum is equal to R„ the symbol rate. However, due to 
the large carrier frequency offset present in the Inmar- 
sat-C signalling and message channels (±1500 Hz), a 
relatively wide bandwidth anti-aliasing filter is required 
at the inputs to the demodulators. For the 1200 bit/sec 
mode of operation (C/N 0 =35.5 dB-Hz), the required 
one-sided filter bandwidth is 2700 Hz, while for the 600 
bit/sec mode (C/N 0 =32.5 dB-Hz) this bandwidth can 
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only be reduced to 2100 Hz. The carrier-to-noise power 
is given by: ftnA 



x Wo b on the unit circle is given by: 

where B is the one-sided filter bandwidth. With the ^"iT^^-lti 1 h*' ' ' 

operating EyN 0 =4.7 dB for both modes and the filter 1* ) 

band widths as above, the carrier-to-noisc power ratios m i_- u * * i * * r •* i rrywn\ 

4 . . . , , . , \ ' \1 ^ ' . , 10 which is equivalent to a finite impulse response (FIR) 

at the demodulator input are 1.2 dB for R,= 1200 bit/- filtcri 0 £ cration . Tf the extra coarse frequency esti- 

sec, and -0.7 dB for R,=600 bit/sec. Therefore, a loss male ^ exacl( this ^ ^ equivalent to a matched filter- 

of approximately 2 dB in the acquisition performance of ing operation for unfiltered BPSK transmission. The 

the 600 bit/sec mode compared to the 1200 bit/sec flher response is given by h(n)=z|-" t 

mode is expected. 15 fj^in^N— 1. 

The frequency acquisition algorithm in accordance With zi-»=e/ 2 ' r * /;v , the frequency response of this 

with the first embodiment of the invention, as has been filter is given by: 
described, consists of two Discrete Fourier Transform 
(DFT) stages. The first stage (i.e., the coarse frequency 

estimator) consists of 64-point hopping Fast Fourier 20 w(<^) « ' c ^ 2irkn/ ^2- n I 

Transforms (FFT) which operate on the squared input n,=0 '*-<A" 

samples. As the carrier to noise ratio (C/N) is reduced 

below 0 dB, the squaring loss increases exponentially. which has the form 

This loss, combined with the 2 dB drop in the input 2J 

C/N due to the wider filter bandwidth, would limit the sin(Mj/2) 

performance of the frequency acquisition algorithm in $in(u/2 + vk/N) 

the 600 bit/sec mode. Accordingly, in accordance with . - - r , . _ _ t 

the second embodiment of the invention, another fre- ^ Wldlh of lhc raam Iobe of lhe res P onsc 15 

quency acquisition algorithm has been devised in order 30 

to improve the demodulator performance at 600 bit/sec. / 2lr ^ s \ 

This embodiment now will be described. <»width = 2 1 "3\T" J 

The second embodiment also uses the DFT/FFT 
concept, as in the first embodiment. In order to reduce ThereforCt the DFT can be represented by k FIR filters, 
the noise bandw.dth of the input to the coarse fre- 35 each haying ^ above centered at kfj/N 
quency estimator, an extra coarse frequency estimator, The inputs to the demodulator (f, =4800 Hz) are 
providing earner frequency estimates to within ±150 windowed by a Hamming window and then are pro- 
Hz, is utilized. Operation of the extra coarse frequency cessed by the above DFTs, with each DFT hopped by 
stage does not require squaring of the samples. This two samples. Hence, the outputs of the DFTs are equiv- 
permits reliable operation for this stage at negative alent to the input signal filtered by bandpass filters of 
values of C/N. After the extra coarse estimation stage, total width 1200 Hz, centered at the DFT frequencies 
the received signal spectrum is processed by a narrow- (i.e., 0 Hz, ±600 Hz, etc.) and then decimated 2:1. The 
band filter of bandwidth 1200 Hz, centered at the esti- hopping and decimation steps are desired so as to re- 
mated carrier frequency. After the filtered data is dc- 45 duce the sampling rate and hence the computational 
modulated by a locally generated carrier at the esti- requirements of the subsequent stages. A minimum out- 
mated frequency, it is input to the currently imple- put sampling rate of 2400 Hz is required to avoid alias- 
mented frequency acquisition stages. The filtering oper- ing. 

ation increases the C/N at the input to the second stage If the signal is centered exactly at one of the DFT 

and hence the loss associated with the squaring opera- SO frequencies, the shape of the filter is matched to the 

tion of the second stage is less savere. The introduction si 8 nal spectrum. In general however, a small mismatch 

of the extra coarse stage (which would precede the will occur which is ofno concern since only the peak of 

coarse stage 220 in FIG. 2) reduces the search range of the spectrum needs to be located, and since no further 

the next stage. Hence, the computational requirements processing on this spectrum will be done. Also, intro- 

should not be worse than in the first embodiment. Fur- » auction of a Hamming window results in improved 

ther details on the algorithm now will be provided. Performance over a rectangular window. The Ham- 

n, * 1 « f . .;. . mine wmdow has a wider transition region and rusher 

Block diagrams of the frequency acquisition schemes . . f . ... . 5 , f 

■ ... . - , . j ... , sidelobe attenuation than the rectangular wmdow. 
,n accordance with the first and second embodiments H in of a mismatch, more of fhe signal spec- 
are shown in FIGS. 7 and 8, respectively. The extra ^ trum yed wUh a Hmmi ^ 
coarse frequency stage of the new algorithm consists of rectangular window 

a series of 8-point DFTs sliding in time. At 600 bit/sec Thc DFT output for the frcquc n C y bin corresponding 

and a sampling rate of 4800 Hz, the frequency bin width l0 lhc maximum magnitude is selected as the likely 

is 1200 Hz. A 300 Hz frequency step size (i.e.. resolu- signal( and then demodulated by a locally generated 

tion) was chosen. Since the frequency range of interest 6 5 carrier with a frequency value corresponding to the 

is ± 1500 Hz, 1 1 outputs need to be calculated. selected bin. Hence, the frequency uncertainty is now 

The spectrum of the N-point sequence, x(n), at point reduced to ± 150 Hz. The subsequent processing stages 

are similar to the ones currently used. The demodulated 
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signal is then squared, windowed by a Blackman win- slightly better than the second embodiment shown in 

dow, and processed through a series of 64-point DFTs. FIG. 8. The trade-off between the two methods also 

Because of the squaring operation, the frequency uncer- involves the faster execution and larger storage space of 

tainty is increased to ±300 Hz. The sampling rate at this the seC ond embodiment, shown in FIG. 8, versus the 

stage is 2400Hz yielding a resolution of 37.5 Hz with a 5 sIower execution and smaller storage space needed for 

64.po.nt DFT. Jwenty^hree DFT outputs are calcu- the „ R flker mc , hod showj] fa pjQ \ (n0 DFT out- 

Ihh , ?° Ver • a ? u " c ^ ,n, y.,r«"8 e ' r h puts need to be stored). Considering the close perfor- 

additional margin. As mentioned before, if the earner f . ^ ' ' 6 " ,. . 

frequency does not coincide with one of the DFT fre- P>ance fthe two echniques at an E^N 0 of 0 dB, which 

quencies of the extra coarse stage, the signal spectrum 10 *"g^*y™<™> the nominal operating E^/N, of 

will be somewhat distorted by the DFTs "filtering" 4/7 dB ' 11 would ap P ear to bc more desirablc to cra P lov 

operation; however, the carrier frequency information the techn,c * ue exemplified in FIG. 8, which employs no 

is still present. external IIR filters. 

The DFT output of the second stage with the maxi- The technique of FIG. 8 also may be implemented 

mum magnitude is frequency corrected by the coarse \$ desirably for the 1200 bit/sec mode in order to improve 

frequency estimate, and input to the fine frequency tne acquisition performance, and minimize the software 

acquisition stage which consists of a single 16-point differences between the two operating modes. With the 

DFT. For this stage the frequency uncertainty is ±37.5 current signalling channel IF test results indicating 

Hz, and the sampling rate is reduced to 75 Hz. The reliable operation at an E*/N<, of 2 dB, introduction of 

frequency resolution for this DFT is about 4.69 Hz, 20 tn ^ new scheme could improve the acquisition perfor- 

Therefore, 17 outputs are calculated. After division by mance of the 1200 bit/sec mode to lower Ej/N 0 values, 

two, to account for the squaring operation, the resolu- The overall structure for the three-stage frequency 

tion is improved to 2.3 Hz. estimator in accordance with the second embodiment of 

The operation of the first two stages of the algorithm the invention is illustrated in FIG. 10. The complex 
were simulated at Ej/N 0 values of 0 and 1 dB. The 25 analog input is filtered by antialiasing filter 1001 and 
results, shown below in Table 2, were obtained by pro- sampled by A/D converter 1002 at a sample rate of F s . 
cessing 500 signalling channel packets with uniformly The estimator requires a block of N complex samples of 
distributed frequency offsets in the range ± 1 500 Hz. a burst transmission stored in its input bufTer 1003. Stage 
The results show that error free frequency estimates are t provides a rough estimate, F, of f f directly from the 
obtained at Es/N„= 1 dB. }Q MPSK-moduIated signal by computing a series of DFT 
TABLE 2 calculations at block 1103, the DFT calculations being 

Simulation Results for the New Frequency Acquisition hopped in time. The time-hopped complex DFT OUt- 

Aigorithm. Using the DFTs "Filtering" Operation puts are bandpass filtered versions of the input signal 

Ej/Nq Noise Seed 100128 Noise Seed 91682 with a two-sided bandwidth of approximately 2 Rs- In 

1 dB Missed o/soo Missed 0/500 35 addition, the time-hopping process decimates the sam- 

0dB Missed 10/300 Missed s/500 pij n g rate to a lower value. Thus the DFT calculation 

simultaneously provides the initial carrier frequency 
Since the misses at an Ej/N 0 of 0 dB were mostly estimate, bandpass filters the input signal, and decimates 
caused by the second (i.e., the coarse) frequency acqui- the sampling rate to approximately 2 R5. The fact that 
sition stage, it seemed worthwhile to determine if this 40 these three operations are combined in this manner 
was in turn caused by the distortion of the received greatly reduces the total number of calculations re- 
spectrum introduced by the first (i.e. the extra coarse) quired. The sample rate decimation also reduces the 
stage. Hence, the simulation program was changed such complexity of the following stages. For an initial uncer- 
that the "filtering'* operation of the first stage was re- tainty f c >R 5 , the bandpass filtering operation improves 
placed by an Infinite Impulse Response (IIR) filter, 45 the signal-to-noise ratio by more than 3 dB. The filtered 
yielding a modification of the second embodiment of samples arc rotated by -F| 1108, reducing the carrier 
the invention. The demodulation operation then was frequency to f f -Fj before being processed in stage 2. 
performed on the delayed input samples instead of the The MPSK modulation is then removed by raising 
DFT outputs. A block diagram of this modified second the stage 1 output samples to the M'th power at block 
embodiment is shown in FIG. 9. The results obtained at 3U 1200. Thus the signal consists of a single carrier at a 
Ej/N o =0 dB for each filter bandwidth are shown in frequency of M(f c -F t ). In a manner nearly identical to 

stage 1, stage 2 provides a higher resolution estimate, 

TABLE 3 F2, of the residual frequency error and simultaneously 

Simulation Results at E^of OdB for the New 55 bandpass filters the carrier signal to a much narrower 

Frequency Acquitition Algorithm. Using the HR Filters bandwidth to permit simultaneous decimation to a 

iir Filter Bandwidth Noise Seed 100128 Noise Seed 916B2 lower sample rate. The samples are rotated by — F2 in 

0.75 kHz Missed 29/500 Missed 18/500 block 1208, reducing the carrier frequency to 

0.5 kHz Missed 8/500 Missed 3/500 f c _F| — fj. 

0 3 kHz Missed 39/500 Missed 40/300 w A high resolution frequency estimate, F 3 , is obtained 

in stage 3 by doing an FFT operation on the stage 2 

Table 3 shows that the acquisition performance de- output samples. Finally the three estimates are summed 

grades if the filter bandwidth is too wide (enough noise at block 1401 to give the final, high -resolution estimate: 
power is not rejected), or too narrow (the signal spec- 
trum is not preserved adequately). 65 

Thus, according to the modified second embodiment j . l ... 

of the invention, with the IIR filter bandwidth of 0.5 F < a Fi + ~ Fj + ~ F > 
kHz, the acquisition performance at an Ej/N 0 of 0 dB is 
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Details of the three frequency estimators are pro- 
vided below. The example packet format of FIG. 1, 
described in detail above, was used for the testing. 

As described above, the frequency acquisition 
method used in the first two stages consists of comput- 
ing a series of time-hopped windowed DFTs (called 
"periodograms"). and averaging the energy magnitudes 
over the duration of a transmission burst to reduce noise 
and fading effects. Given a sequence of N buffered 
complex samples, x(n), where n=0, 1, 2, . . . , N - 1, in 
block 1003, it is possible to form K subsequences x^n), 
r= 1, 2, . . . , K, of length L spaced J samples apart. 

For each of the K subsequences, the windowed DFT 
X^k) is given by the following equation: 



10 



15 



XAk) . 



■j(lit/L)nk 



(2) 



where 

W(n) is the window function 

k = DFT frequency bin number. 
The foregoing equation is computed in block 1103 for 
stage 1 using Kaiser window function 1102. The calcu- 
lation is performed for all values of k within the desired 
frequency range and repeated for all K i subsequences. 

The periodogram \ASk) is equal to the squared magni- 
tude of a DFT bin: 
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where 



DFT frequency. 

The unnormalizcd power spectral density is then: 
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In stage 1 the periodogram I^fjt) is computed in block 
1104 for all Ki subsequences and for all frequency bin 
numbers, k, within the frequency range of interest. 
S*(f*) is computed in block 1105. 

The power spectral density is simply a sum of the 45 
periodograms of each subsequence and the normalized 
carrier frequency is given by f* for the value of k which 
maximizes S^f*). Thus the frequency estimate is ob- 
tained by searching for the maximum S*(f*) for all val- 
ues of k (block 1106 in stage 1). The step size and range 
of k is arbitrary and depends on the application. If the 
DFT is implemented as a Fast Fourier Transform 
(FFT), then k takes on the values 0, 1 L- 1. 

For the example case, N = 2048 samples (256 sym- 
bols) and F;=4800 Hz (8 samples per symbol). Only the 
symbols in the middle of a TDM A timestot assumed to 
contain the carrier arc processed. The remaining sym- 
bols at the beginning and end of the timeslot are dis- 
carded due to burst location uncertainty. For R,=600 
bit/s, the width of the main spectral lobe is 600 Hz. For 60 
a sample rate of 4800 Hz, L» = 8 provides a DFT bin 
width that matches the main spectral lobe of the BPSK- 
modulated signal (FIG. 11). A Kaiser window (block 
1102) with parameter >3 = 3. 1 was chosen for the win- 
dow function to maximize the average signal-to-noise 
ratio within a DFT bin for all expected values of fre- 
quency offsets. The frequency resolution of the first 
estimator is 300 Hz— i.e. k has a step size of 0.5, DFT 
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outputs are computed for all values of k within the 
range of ±5.0 to span the total specified frequency 
uncertainty range of ±1500 Hz. The first stage esti- 
mate, Fj, is kmox-Ff/Ji, where k max corresponds to the 
DFT bin with the maximum power spectral density and 
Ji represents the sample rate decimation ratio of stage 1. 

Each complex DFT output X,(k) can be thought of as 
an output at time r of a Finite Impulse Response (FIR) 
bandpass filter with impulse response W centered at 
frequency k. Thus the sequence Xr(kmax), r= 1, 2, . . . , 
K] is a bandpass filtered version of the input signal and 
the signal-to-noise ratio is greatly improved. Since each 
subsequence is spaced J] input samples apart, the se- 
quence Xf(lL m ax) has a sample rate decimated by a factor 
of Jj compared to the input sequence. XJk max ) is buff- 
ered in block 1107 for subsequent processing by stage 2. 

J I must be chosen carefully to avoid aliasing effects. 
Fj/J| must be greater than the filtered signal bandwidth. 
The two-sided signal bandwidth of 6PSK assuming no 
transmit filtering is approximately 4R,. Therefore Jj = 2 
and the decimated sample rate = 2400 Hz. The sequence 
X^kmox) stored in block 1107 is frequency corrected by 
~F\ to greatly reduce the frequency offset (block 
1108). Frequency correction of the stored complex 
samples can be performed by multiplying the samples 
by [cos(27rFin/F f ) — Jsin(27rFm/Fj)] ( where n is the 
sample index. The sine and cosine values are obtained 
from a look-up table. Precautions need to be taken in 
minimizing the truncation error while calculating the 
argument 27rF(n/F J . 

Next, the modulation is removed in block 1200 by 
squaring the complex inputs (for BPSK modulation). 
This concentrates most of the signal energy at a single 
frequency equal to twice the modulated carrier fre- 
quency. Although the squaring of the signal results in a 
loss in the operating signal-to-noise ratio, this loss is 
regained in the next stage due to the averaging process 
performed on the DFT magnitudes. 

The 1024 complex samples from stage 1 are divided 
into K2=31 blocks of L2— 64 symbols spaced 
^2 = ^2/2 = 32 symbols apart (block 1201). A spacing of 
L/2 produces the most reliable estimate versus the 
amount of computations required. The periodograms 
are computed in a similar manner to stage 1 using equa- 
tions 2, 3, and 4 in blocks 1203, 1204, and 1205 respec- 
tively. Each subsequence is multiplied by a 64 point 
Blackman window in block 1202. A Blackman window 
has a very broad transition width so that a carrier lo- 
cated at the very edge of a DFT frequency bin will be 
attenuated only 1.1 dB as compared to 3.9 dB for a 
rectangular window. At a sampling frequency of 
Fi=2400 samples/sec a 64 point DFT provides a reso- 
lution of 37.5 Hz. The resolution is increased to 18.75 
Hz since the frequency estimate of the squared data is 
divided by two in block 11209. (In general, removing the 
modulation of M-ary PSK by raising the data samples to 
the Mth power increases the carrier frequency by a 
factor of M.) Because the stage 1 frequency resolution is 
300 Hz, the stage 2 search can be limited to ±150 Hz. 
However, the search is extended to ±412.5 Hz (DFT 
frequency bins in the range k= ± 1 1) to compensate for 
estimation error in stage 1. A DFT can be implemented 
more efficiently over the limited search range than the 
FFT implementation. The frequency bin with the high- 
est energy, kmu. is assumed to contain the carrier and 
the estimate, F2, is determined (block 1206). The stage 2 
sequence X^k mox ) can also be viewed as a bandpass 
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filtered and decimated version of the stage 2 input data 
and is buffered in block 1207 for further processing. The 
output sample rate is decimated by a factor of L2/2 to 
75 Hz. The 32 samples stored in block 1207 are fre- 
quency corrected by -F 2 in block 1208 and then are 5 
further processed by stage 3. 

The stage 2 output samples are truncated to one block 
of L3 = 16 complex samples (block 1301) and a 16-point 
FFT calculation is performed in block 11302. The FFT is 
more efficient than a DFT in this case. The squared 10 

magnitude of all 16 FFT outputs (k=0, 1 15) is 

calculated in block 1303. The frequency bin with the 
maximum magnitude is selected in block 1304 and the 
frequency estimate F3 is determined with a resolution of 
4.69 Hz. The actual resolution is 2.34 Hz since the esti- 15 
mate is divided by two (block 1305) to compensate for 
the squaring operation. For mobile satellite applica- 
tions, the actual resolution is limited by short-term Dop- 
pler effects caused by motion of the mobile terminal and 
is on the order of 5 to 10 Hz. 

The total frequency estimate then is the sum of the 
three individual estimates, as given in equation (1). Note 
that each stage attempts to correct for frequency errors 
in the previous stage by using an extended frequency 
search range. 

The frequency estimator will always produce an out- 25 
put regardless of the presence of a carrier. The carrier is 
assumed to reside in the DFT frequency bin with the 
highest average energy. A good measure of the reliabil- 
ity of the estimate can be easily determined by the peak- 
/average energy ratio, R, from stage 1 or 2. Normally a 
strong spectral peak is observed in one frequency bin 
and low-level noise is observed in all the other bins. The 
peak energy is simply the average energy of the peri- 
odogram containing the maximum energy (Sj((f* m „)). 
The average energy is the average of all the periodo- 35 
grams for all other values of k. However, since signifi- 
cant carrier energy leakage can be observed in fre- 
quency bins adjacent to k maJt , adjacent bins also should 
be excluded from the average. Thus: 
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SAfkmax) 



kz£k max ,k-#.k mog — l,k^.k max + I 



This ratio is also a good means of detecting transmis- 45 
sion burst activity. If the ratio is greater than a predeter- 
mined threshold, then a burst has been detected. This 
method is very reliable over a wide dynamic range, 
eliminating the need for gain control. Because a ratio is 
computed, a receiver gain factor would appear in both 50 
the numerator and denominator and therefore be can- 
celled out. Simulations indicate that a peak/average 
energy ratio threshold near 2.0 produces good results 
down to Ei/N 0 =2 dB. 

While the invention has been described above in de- 55 
tail, various changes and modifications within the scope 
and spirit of the invention will be apparent to those of 
working skill in this technological field. Thus, the in- 
vention is to be considered as limited only by the scope 
of the appended claims. 60 

What is claimed is: 

1. A digital preamblclcss demodulator comprising: 
an input buffer for receiving an input signal which 
includes a carrier wave and providing a baseband 
output; 65 
a coarse frequency estimator, connected to said input 
buffer, for providing a first estimation of a fre- 
quency of said carrier wave, said first estimation 



being within a first frequency band, said coarse 
frequency estimator comprising a first Fourier 
transform module for performing either a first fast 
Fourier transform (FFT) or a first discrete Fourier 
transform (DFT) to provide said first estimation; 
and 

a fine frequency estimator, receiving said first estima- 
tion and providing a second estimation of said fre- 
quency of said carrier wave, said second estimation 
being within a second frequency band that is nar- 
rower than said first frequency band, said fine fre- 
quency estimator comprising a second Fourier 
trans form module for performing either a second 
FFT or a second DFT on said first estimation to 
provide said second estimation; 

the foregoing arrangement being such that said fre- 
quency of said carrier wave can be acquired within 
a predetermined resolution without necessitating a 
preamble in said input signal. 

2. A digital preambleless demodulator as claimed in 
claim 1, wherein said first frequency band is ±75 Hz. 

3. A digital preambleless demodulator as claimed in 
claim 1, wherein said second frequency band is ± 10 Hz. 

4. A digital preambleless demodulator as claimed in 
claim 1, further comprising an extra coarse frequency 
estimator, disposed between said coarse frequency esti- 
mator and said input buffer, for receiving said baseband 
output and providing a third estimation of said fre- 
quency of said carrier wave, said third estimation being 
within a third frequency band that is wider than said 
first frequency band, said extra coarse frequency estima- 
tor comprising a third Fourier transform module for 
performing either a third FFT or a third DFT on said 
baseband output to provide said third estimation, said 
coarse frequency estimator receiving said third estima- 
tion and providing said first estimation accordingly. 

5. A digital preambleless demodulator as claimed in 
claim 4, wherein said third frequency band is ± 150 Hz. 

6. A digital preambleless demodulator as claimed in 
claim 5, wherein said first frequency band is ±37.5 Hz. 

7. A digital preambleless demodulator as claimed in 
claim 5, wherein said second frequency band is ±2 3 
Hz. 

8. A digital preambleless demodulator comprising: 
an input buffer for receiving an input signal which 

includes a carrier wave and providing a baseband 
output; 

an extra coarse frequency estimator for receiving said 
baseband output and providing a first estimation of 
a frequency of said carrier wave, said first estima- 
tion being within a first frequency band, said extra 
coarse frequency estimator comprising a first Fou- 
rier transform module for receiving said baseband 
output and performing a first fast Fourier trans- 
form (FFT) or a first discrete Fourier transform 
(DFT) on said baseband output to provide said first 
estimation; 

a coarse frequency estimator, connected to said extra 
coarse frequency estimator, for providing a second 
estimation of a frequency of said carrier wave, said 
second estimation being within a second frequency 
band that is narrower than said first frequency 
band, said coarse frequency estimator comprising a 
second Fourier transform module for performing 
either a second FFT or a second DFT on said first 
estimation to provide said second estimation; and 
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a fine frequency estimator, receiving said second input signal to provide a first estimation of said 

estimation and providing a third estimation of said frequency of said carrier wave, said first estimation 

frequency of said carrier wave, said third estima- being within a first frequency band; and 

tion being within a third frequency band that is performing a second FFT or a second DFT on said 

narrower than said second frequency band, said 5 first estimation to provide a second estimation of 

fine frequency estimator comprising a third Fou- frequency of said carrier wave, said second 

nertransform module for performing either a third estimation being within a second frequency band 

l!7 ■? ' ^".i a . ° n SCC ° nd es,imation 10 that is narrower than said first frequency band. 

,hf CI" , " tm «™; fc u . A , , n 13. A method as claimed in claim 12, wherein said 

JZ,I X B a I ran * emen, bem « ^ ucn ,ha, ^ ud l re - 10 first frequency band is ±75 Hz. 

T^^l^iS^^^, 7 h ' n »*• A "nethod as claimed in claim 12, wherein said 

a predetermined resolution without necessitating a 

preamble in said input signal. g ~ d A fre ™ b ? d 10 , H \, , . 

9. A digital preambleless demodulator as claimed in . I5 ' A mctho f planned m claim JZj furthei ' compns- 
claim 8, wherein said first frequency band is ± 1 50 Hz. 1 < ,ng the steps of Panning a third FFT or a third DFT 

10. A digital preambleless demodulator as claimed in on M,d m P ut SI « naI 10 P rovide a third estimation of said 
claim 8, wherein said second frequency band is ±37.5 frequency of said earner wave, said third estimation 
Hz. being within a third frequency band which is wider than 

11. A digital preambleless demodulator as claimed in first frequency band, said first estimation being 
claim 8, wherein said third frequency band is ±2.3 Hz. 20 provided accordingly. 

12. A method of performing digital preambleless 1& A method as claimed in claim 15, wherein said 
demodulation of an input signal, said method cornpris- first frequency band is ±150 Hz. 

ing: 17. A method as claimed in claim 15, wherein said 

receiving an input signal which includes a carrier second frequency band is ±37.5 Hz. 

wave; 25 18. A method as claimed in claim 15, wherein said 

performing a first fast Fourier transform (FFT) or a third frequency band is ±2.3 Hz. 

first discrete Fourier transform (DFT) on said * * * • • 
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